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Abstract 

We present some exact results about universal quantities derived from the local 
density matrix p, for a free massive Dirac field in two dimensions. We first find tvp n 
in a novel fashion, which involves the correlators of suitable operators in the sine- 
Gordon model. These, in turn, can be written exactly in terms of the solutions of 
non-linear differential equations of the Painleve V type. Equipped with the previous 
results, we find the leading terms for the entanglement entropy, both for short and long 
distances, and showing that in the intermediate regime it can be expanded in a series 
of multiple integrals. The previous results have been checked by direct numerical 
calculations on the lattice, finding perfect agreement. Finally, we comment on a 
possible generalization of the entanglement entropy c-theorem to the alpha-entropies. 

1 Introduction 

The trace of the vacuum state projector over the degrees of freedom corresponding to a 
spatial region, results in a mixed density matrix with a non-vanishing 'geometric' entropy. 
This kind of construction was proposed by some authors ^ |2J |3J 11] in an attempt to 
explain the black hole entropy as some sort of entanglement entropy of the corresponding 
vacuum state. However, since the role of gravity cannot be ruled out, this identification 
is still a conjecture. 

In recent years, there has been a renewed interest on the properties of local reduced 
density matrices, specially for low dimensional systems (see for example E3 El El EH EH 
I11U12| . and references therein). This was partially motivated by developments in quantum 
information theory and also on the density matrix renormalization group method in two 
dimensions ^3] • These investigations have made it manifest that the entanglement entropy, 
as well as other measures of information for the reduced density matrices of the vacuum 
state are interesting quantities by their onw right. Besides, they can yield a different view 
on certain aspects of quantum field theory, since the objects one has to calculate are quite 
different to the standard ones. That means, for example, that there is a different structure 
of divergences, and a nice interplay between geometry and UV behaviour. 

In this paper, we study measures of information for the local density matrix pa for a 
Dirac fermion in two dimensions. This is achieved by tracing the vacuum state over the 
degrees of freedom outside the set A. Specifically, we shall consider the a-entropies 

S a (A) = _L- logtr(^), (1) 
1 — a 
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and the entanglement entropy 



8(A) 



tr(p A log pa) = lim S a (A). 



(2) 



These two functions of pa are perfectly well defined on a lattice, but their continuum 
limit are plagued by UV divergences. In two dimensions, however, the form of those 
divergences is particularly simple, being just an additive constant proportional to the 
logarithm of the cutoff and to the number of boundary points in A. Indeed, since the 
divergences have their origin in the UV fixed-point (where masses can be neglected), this 
follows from the result for the conformal case E] • Remarkably, for the entanglement 
entropy, it also follows from the strong subadditive property of the entropy, when the 
spatial symmetries are taken into account 6 . 

Several universal quantities can be easily obtained from S(A) and S a (A). In particular, 
when A is a single interval of length r, we define the dimensionless functions 



The function c(r) is always positive and decreasing, and we shall call it 'entropic c- 
function', since it plays the role of Zamolodchikov's c-function in the entanglement entropy 
c-theorem j^j. All the universal information which can be obtained from S a (r) and S(r) 
is encoded in c a (r) and c(r). 

For more general sets, formed by several disjoint intervals, universal quantities can be 
constructed through the mutual information function /, which for two non intersecting 
sets A, B is given by 



Remarkably, / remains finite in higher dimensions 1 . 

The traces of powers of p involved in the a-entropies, with a = n £ Z, can be repre- 
sented by a functional integral with the fields defined on an n-sheeted surface with conical 
singularities located at the boundary points of the set A. This integral is quite diffi- 
cult to deal with, except for the conformal case, where that surface can be conveniently 
transformed. 

We shall use here a novel approach (see section |2J) , whereby the problem in the n- 
covered plane is mapped to an equivalent one in which an external gauge field couples to 
an n-component fermion field (defined on the plane) and its role is to impose the correct 
boundary conditions through its vortex-like singularities. Then the resulting theory is 
bosonized, to express trp n as a sum of correlators of local operators. 

In the massless case, discussed in section El the bosonized theory becomes a massless 
scalar field, and thus we could obtain the entropies explicitly, with results in agreement 
with the ones of 0. As a by-product we derived a quite convenient expression for I (A, B). 

For the massive case (section |1J) , the dual theory is instead a sine-Gordon model at the 
free fermion point. In this case, we still could evaluate trp n exactly for a single interval 
of size r, by expressing it as a sum of n correlators for exponential operators. A method 
to find that type of correlators was (fortunately) already available in the literature |15j . 
and it could be adapted to our case after some minor modifications. The outcome is an 

1 Formulae Q make sense only when a translation invariant cutoff is used. This appar- 
ent limitation could be overcome, in two dimensions, by defining c(r) through the function 
F(A, B) = S(A) + S{B) - S(A n B) - S{A U B) used in 6 for overlapping sets. However, F(A, B) is 
not finite, in general, for dimensions greater than 2. 




(3) 




(4) 
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Figure 1: The plane with cuts along the intervals (ui,Vi), i = l,...,p. The circuits C Ui and 
C Vi are used in the text to discuss the boundary conditions. 

exact expression for trp n in terms of the solutions of second order non linear differential 
equations. Also for the massive fermion case, we obtained expansions for small and large 
values of mr (m is the fermion mass), and by analytical continuation in n we calculated 
the corresponding expansions for the entropy. 

This article concludes with a summary of our results (section |SJ) and a discussion on 
a possible generalization of the c-theorem to the a-entropies. This conjecture naturally 
suggests itself by the results of the explicit calculations. 



2 Reduced density matrix for a Dirac fermion 

We consider the reduced density matrix for a free Dirac fermion, which is defined by 
tracing the vacuum state over the degrees of freedom lying outside a given set A. We 
specify A as a collection of disjoint intervals (ui,Vi), i = 1, ...,p (see figure 1). 

Following |14j . the density matrix p{^m, ^out) can be written as a functional integral 
on the Euclidean plane with boundary conditions \P = \Pi n , and \P = ^out along each side 
of the cuts (ui, Vi) 

p(*in,*out) = ^TT / D*e- s W, (5) 



Z{1] 

where Z[l] is a normalization factor, introduced in order to have trp = 1. 

In order to obtain trp n , we consider n copies of the cut plane, sewing together the cut 
(iii,t>i)o Ut w ith the cut (uj, v i)f n +1 , for all i = 1, ...,p, and the copies k = 1, ...,n, where the 
copy n + 1 coincides with the first one. The trace of p n is then given by the functional 
integral Z[n] for the field in this manifold, 

For a fermion field, we have to remember that the trace requires to introduce a minus 
sign in the path integral boundary condition between the fields along the first and the 
last cut ^H]) as it happens for the fermion thermal partition function in the Matsubara 
formalism ^7j. Besides, in the present case, for each copy there is an additional factor 
— 1. This is due to the existence of a non trivial Lorentz rotation around the points Ui and 
Vi which is present in the Euclidean Hamiltonian when expressing trp n as a path integral 
|18| . From these considerations, we finally get a total factor (— l)( n+1 ) connecting the 
fields along the first cut (ui,Vi)j n and the last one (^i,Wj)out i n @- 

Rather than dealing with field defined on a non trivial manifold, we find it more 
convenient to work on a single plane, although with an n-component field 



/ *x(x) 



(7) 
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where ^i(x) is the field on the I th copy. Of course, the singularities at the boundaries are 
still there, and we shall show a simple way of taking them into account now. 

Note that the space is simply connected but the vector ^ is not singled valued. In fact, 
turning around any of the C Ui circuits (see figure 1) it is multiplied by a matrix T, and 
after turning around the C v . circuit it gets multiplied by the inverse matrix T _1 . Here, 



/ 1 \ 

1 

1 

V J 



(8) 



which has eigenvalues e*« 27r , with k = — ^ 2 ^ , — ^ + l,... / n 2 1 ^ . 

Then, changing basis by a unitary transformation in the replica space, we can diag- 
onalize T, and the problem is reduced to n decoupled fields <fr k living on a single plane. 

■ k o 

These fields are multivalued, since when encircling C Ui or C Vi they are multiplied by e % ~ 77 

■ k r> 

or e~ l ™ 7r , respectively. 

That multivaluedness can now be easily disposed of, at the expense of coupling singled- 
valued fields 3> fc to an external gauge field which is a pure gauge everywhere, except at 
the points Uj and Vi where it is vortex-like. Thus we arrived to the Lagrangian density 

C k = !> V {df, + i A k ^j $ k + ml> fc $ fc . (9) 

Indeed, the reverse step would be to get rid of the gauge field by performing a singular 
gauge transformation 

$*( x ) _> e -*/"o dx'"Ai(x')^ k (x) ; (1Q) 

(where xq is an arbitrary fixed point). Since the transformation is singular, one goes 
back to a multivalued field. By the same token, we learn that in order to reproduce the 
boundary conditions on <£ fc , we should have 

flr'-AtAx) = -—, (11) 
n 



\k 



2irk 



dx»A*Jx) = . (12) 



n 



As a side remark, we note that the addition of terms of the form 2-7rg, with q an integer, to 
the right hand side of the above formulae does not change the total phase factor along the 
circuits. However, the winding number of the phase, and thus the boundary conditions 
the gauge field imposes, would be different. In fact, the free energy does depend on these 
integers and choosing q ^ would not select the vacuum state. 

Equations (|11[) and (|12|) hold for any two circuits C Ui and C v . containing Uj and Vi 
respectively. Thus: 

e^d v A k (x) = 2^ K x " u i) ~ 6 ( x ~ v $ > ( 13 ) 
i=l 

where the presence of a vortex-antivortex pair for each k and each interval is explicit. The 
functional integral is then f actor ized: 

(n-l)/2 

Z[n]= [] Z k , (14) 
fc=-(n-l)/2 
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where Z k can be obtained as vacuum expectation values in the free Dirac theory 

Z k = (e^ A ^\, (15) 



where j£ is the Dirac current, A 1 ^ satisfies ()13p . and we adopted a normalization such that 
<1> = 1. 

3 Bosonization and the massless case 

In order to evaluate ()15|) . it is quite convenient to use the bosonization technique |19| . to 
express the current j£ as 

j^-L e ^d v <f>, (16) 



where is a real scalar field. For a free massless Dirac field, the theory for the dual field 
(j) is simply 

C^^d^cp. (17) 

Therefore we have to evaluate: 

Zk = / e ijA%j s d»'9 v <tf'x\ = / e -iV5FiES=iW«i)-*(«i))\ , (18) 



where the vacuum expectation values correspond to the theory (|17j). Since is quadratic 

1 ' e -iJ f(x)<K*)#*\ = e -|/ f(x)G(x-y)f(y)d^xd 2 y ^ ^ 



with the correlator 

G(x-y) = -^-log\x-y\ , (20) 
it follows that (|18|) can be written as 

2k 2 

\ogZ k = T E(ui,Vj) , (21) 

E(ui,vj) = y~] log \uj -Vj\—y~] log k - uj\ - y~] log - - piogg . (22) 

Here e is a cutoff introduced to split the coincidence points, \ui — Ui\, \vi — Vi\ — > e. 
Summing over k and using Q and (j2J we obtain 

S„ = log(trp») = — ?— ^logZ fe = i^E(^), (23) 

1 — n 1 — n b n 



5 = -S (24) 



This agrees exactly with the general formula for the entanglement entropy for conformal 
theories obtained in [3] . The general case differs from Q24[) on a global factor of the Virasoro 
central charge C, where C = 1 for the Dirac field. 

Equation (|24|) has an interesting corollary: recalling the definition Q for the mutual 
information, it follows that, for non-intersecting sets A, B and C 

I(A,BuC) = I(A,B) + I{A,C). (25) 
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That is, in contrast to the entropy, the mutual information is extensive (in each of the 
sets separately) in the conformal case. This curious property does not hold in the non 
confer mal case or in more dimensions. It can be written as 



3/4 Ib y {x - yY 



I(A, B) = - I dx I dy- ? (26) 



for any two non intersecting A and B. This formula for I(A,B) shows explicitly its 
model independent properties, that is, it is cutoff independent, positive and monotonically 
increasing with A and B. 

For a single interval we have that 

a + 1 1 
c a = — — , c = -, (27) 
ba 6 

which are constants. 



4 The massive case 

We consider here a massive fermion, focusing on the case of a single interval of length r. 

We can still use (|16|) to deal with l|15|) for the partition function. Now, however, the 
bosonization of the massive fermion theory leads to a sine-Gordon Lagrangian ^S] 

C = ^ (d^cj) + A cos(\/I^>)) , (28) 

where A is a mass parameter. Then, as in the previous section we have 

(n-l)/2 

log(trp n (r)) = logU'^^e^^) , (29) 

fc=-(n-l)/2 

where now the expectation value is evaluated in the sine-Gordon theory at the free fermion 
point given by the Lagrangian l|28jl . The correlators of exponential operators were studied 
in ref. [T5|, where it was shown that | e 'V5<»(t(f) e >v / iaV(0)j f or a ^ a ' G [ , 1] can be 
parametrized by a function satisfying a nonlinear second order differential equation of the 
Painleve type. This is done by mapping the sum over the form factors into the determinant 
of a Fredholm operator. The general relation of this type of determinants and differential 
equations is studied in [20] . The result of ^3120] can n °t be directly applied here, since we 
need the correlator for a = —a. This case can be dealt with through a minor modification 
in those results (see Appendix A for more details). 
To proceed, one introduces the function 

Wa ( x ) = r j- log ^ e -V4^(r) e iVW(0)^ ^ (3Q) 



with 



00 

2 



w a (x) = - yv a (y)dy, (31) 



where v a satisfies 

x 1 — K \ / x z 1 — v a z 
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Figure 2: Solid lines are plots of c n (x), obtained by solving differential equations Q31I34JI 
numerically. The values of n are, from top to bottom, n = 2, 3, 5 and 50. These functions 
take the value (n + l)/(6n) at the origin, which cumulates at 1/6 = 0.166... for large n. For 
large x, they decay exponentially fast. Dotted lines correspond to the c n (x) that results 
from putting the model on a lattice, for n = 2 and n = 3. This points are evaluated for set 
sizes ranging from 200 to 600 lattice points and inverse mass values ranging from 200 to 
3200 lattice units. The fact that these points, computed for several different lattice mass 
values, tend to lie in a single continuous curve shows the universal character of c n (x). 

Here we have defined x = rm, and the boundary condition for H32|) is 

2 

v a (x) ~ — sin(7ra)K2a(^) as x — ► oo , (33) 

7T 

where K2 a (x) is the standard modified Bessel function. 

Thus, (|3*Tj) . (|3*2*)l and ()33|) give the exact value of c n (r) in terms of a finite number of 
Painleve type functions 

(n-l)/2 

CnO) = z V w k / n (x) . (34) 

1 — n ' 
fc=-(n-l)/2 

The functions c n (x) are shown in figure (2) for some values of n. They take the 
massless case value at x = 0, and lie between 1/6 for n = oo and 1/4 for n = 2. At 
large x = rm they decay exponentially fast. Note that the case n = 1 corresponds to the 
entropic c-function, which at the origin takes the value 1/3 (see figure 3). 

We have also made a direct numerical evaluations of c n (x) on the lattice, with a method 
which is described in Appendix B. Figure 2 shows the results for n = 2 and n = 3, which 
match the exact theoretical values given by (|34jl. For higher values of n, the approach to 
the continuum limit of c n {r) is slower, and agreement within few percent requires lattice 
sets bigger than 1000 points, already for n = 4. 

In order to evaluate c n (r) for non integer n, or to compute the entropic c function, 
we would need an analytical continuation in n of (|34|) . In what follows we obtain these 
quantities for expansions at long and short distances. 
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4.1 The long distance expansion 

A naive expansion for large x in (|33j) is not sufficient to evaluate the entropy, which requires 
the a — * 1 limit of c a . Indeed, expansions for large x can be obtained from (|33|) by using 
the asymptotic form of the Bessel function for large values of its argument 



K a (x) 



This gives for the sum of v 2 



(n-l)/2 

2 '■/.-,■ 

fc=-(n-l)/2 



7T 1 



+ 



7Tl 



(4a 2 -1)— + 



.r 2 



2\ 2 (n ^ /2 fc 

-) V sin 2 (^-)K 2 fc (r) 

7 fc=-(n-l)/2 



-2,c 



1 



rn 1 

— + I™ — 

tt x 12 irn 



4 + 2cos(f) 



7rnsin 2 (r)' £ 



- + 



(35) 



(36) 



Now we can apply this formula for non integer values of n = a. This provides a good 
expansion for fixed a and sufficiently large x, from which the expansion for c a (x) follows 
by term by term integration 



-2x 



c a {r) 



a 



a , 1 a 2 — 4 
2^ + ^24 vra 4 



7rasin 2 (-) x 



(37) 



However, we see from the way the coefficients in the expansion depend on a that the series 
does not converge unless x 3> (a — l)" 1 for a approaching 1. 

Thus, (|37jl cannot be used to compute the entropy. This can be repaired using di- 
rectly the first term in the form factor expansion for the correlators (see Appendix A) or, 
equivalently, the integral representation for the Bessel function 



K a (x) 



du e 



+ vV - 1 ) +(u + V 



1 



2v^? 



1 



(38) 



in expression (|36|) and summing over k inside the double integral. We have the following 
asymptotic expression 

e -y(u+v) 



c a (x) 



~~ 9 / dyy 



7T- 



x F n 



du 



dv 



\Ju 2 



1 



(v + Vv 2 -l))+F a ( 



+ vV - 1 



where 



1 



4(1 - a) 




2 cos (|) 



2« 



2_ _ 2 



2cos(f) 



+ 



_ 1_ 

Z a 



(39) 



Surprisingly enough, this function is proportional to a delta function in the limit a. — * 1. 
Indeed, for a fixed £ ^ 1, F a (z) — ► as a — ► 1. However, for z = 1 this limit is singular. 
Around z = 1 and a = 1 the singularity can be isolated by using polar coordinates in the 
(a — 1, z — 1) plane and expanding in the radial coordinate. The result is 



FJz) 



Tr- 



ee 



1) 



2 vr 2 (a - l) 2 + (z - If 



(l + 0((2-l),(a-l))) . 



(40) 



S 




Figure 3: The dotted curve is the function c(r) evaluated on a lattice, with points obtained 
for different values of the mass and lattice distance (see caption of figure (2)). The conti- 
nuity of the plot agrees with the universal character of c(r), which is in fact a function of 
x = mr. The solid- line curves are the short and long distance leading terms we evaluated 
analytically. 

Thus, we have 

lim F a (z) = ~S(z - 1) , (41) 

a— »1 Z 

which yields 

r°° i 

c(x) ~ - J dyyK (2y) = -xK 1 (2x), (42) 

in perfect agreement with the numerical results (see figure (3)). 

This term corresponds to the first one in the form factor expansion, which is due to 
one soliton-antisoliton pair. In the fermion language it is the contribution from a single 
fermion-antifermion pair. It is interesting to note that (|42j) is 1/4 at the origin, and thus 
a single fermion-antifermion pair contributes more than 75% to the entanglement entropy 
function c(r) at all scales. 

More terms in the expansion for the entropy can be written in terms of multiple 
integrals, by using the form factor series for the correlators (|55|). First we expand the 
logarithm of the correlator order by order, and then sum over k inside the integrals. This 
is easily done since we have a formula analogous to (|41j) 

lim > sin 29 -7T )z- = . 2J 5 (z - 1 . 43 

n^il-n ^ \n J 2T(q) y ' K ' 

fc=-(n-l)/2 v 7 yHI 

We do not write explicitly more terms of the series here, since the expression for the multi- 
ple integrals, though straightforward to obtain, do not seem to be particularly illuminating. 
Besides they are difficult to evaluate numerically. 
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4.2 The short distance expansion 



Close to the conformal limit, the best way to expand v a (x) is by a direct use of the 
differential equations. We have the series solution of (|32|) around the origin 

v a (x) = -2alog(x) + 6 + x 2 (2a - 8a 3 + b - 8a 2 b - 4ab 2 - b 3 ) + 

+\ (-a + 8a 3 + 8a 2 b + 3ab 2 ) log(x) 

+ (-4a 3 - 3a 2 o) log 2 (x) + 2a 3 log 3 (x)) + 0(x 4 ) (44) 
which is of the general form 

oo 2s+l 

v a {x) = ^x 2B ^f a , t log t (x). (45) 

s=0 t=0 

The full expansion requires only the knowledge of the constant term /o,o = b. It might be 
absorbed into the logarithm, replacing log(x) by log(xe~2S') and setting 6 = everywhere 
in the above expression. Unfortunately we do not know the general expression for b as a 
function of a. An exception is the case a = 1/2, where b = — 1 — ^e + 31og(2), and ^e 
is the Euler constant. This follows from the results for the Ising model correlators [2TJ 
(see Appendix A). Here we consider only the leading term which is independent of b. The 
integration constant for w a is given by the conformal limit of section 2, w a (0) = —2a 2 . 
Thus we have 

Wa = -2a 2 + 2aVlog 2 (x) + C(x 2 log(x)) , (46) 

c a (x) = ^±i(l-x 2 log 2 (x))+0(x 2 log(x)), (47) 
oa 

c(x) = ^-^x 2 log 2 (x) + 0(x 2 log(x)). (48) 



5 Final remarks 

To summarize the results, we have found the exact function c a (r) for a massive fermion 
field in two dimensions for integer values of a = n. It is given in terms of the solutions of 
non linear differential equations. This completely determines the a-entropies for integer a, 
except for a non-universal, ultraviolet divergent, (additive) constant. We have also found 
the leading terms of the entanglement entropy c-function c(r) for short and large distances 
showing how it can be expressed as a series of multiple integrals. In the massless case, we 
have calculated the entanglement and the a-entropies exactly for an arbitrary set, with 
results that coincide with the ones of 0- Surprisingly, the mutual information turns out 
to be extensive in the massless case. 

There are also other interesting universal quantities which can be derived from the 
ones above. For example, the a-entropies increase logarithmically in the conformal limit, 

S a (r) = ^log(r) + k , (49) 
ba 

and saturate for rm 1, 

S a (r) (50) 
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Although the saturation constant fcoo is cutoff dependent, its dependence on the mass is 
universal and can be expressed in terms of c a . Defining an interpolating function 

k(r) = s a (r) - Ca(r) log (r) , (51) 

we have k(0) = ko, fe(oo) = k^. Using k'(r) = — d a {r)\og{r) it follows that 

/■OO pQC 

Ak = / k'(r)dr = — / log(r)c^ (mr)d(mr) = —Ac a log(m) + const. . (52) 
Jo Jo 

This is the result obtained in from the analysis of the properties of the energy momen- 
tum tensor on a conical space. 

There is a very interesting point arising from the equation 1)3 ljl which gives w a as an 
integral of an explicitely negative quantity. This implies through equation that the 
dimensionless functions c n (r) are always positive and decreasing with r for any integer 
n > 2. On the other hand, the entanglement entropy c-theorem shows that c(r) = c\(r) is 
a dimensionless, positive, and decreasing function for any two dimensional theory [H]. This 
seems to suggest that the same is to be expected for the alpha entropy c- functions c a (r) 
for any a. This is equivalent to say that S a (ri) — S a (r 2 ) is decreasing under dilatations 

for n > V 2 , 

S a (\ri) - S a (Xr 2 ) < S^n) - S a (r 2 ) , (53) 
with A > 1. The corresponding relation for the entanglement entropy (a = 1) follows from 

El 

Remarkably, in that case each c a would lead to an alternative form of the c-theorem 
in two dimensions. Note that the c a (r) are universal quantities which have fixed point 
values proportional to the Virasoro central charge. 

Recently, a different conjecture which also involves the renormalization group flow and 
the reduced density matrices has been presented in a series of papers [313101 HQ] . Those 
authors proposed several interesting majorization relations for the local density matrices 
and connect them to the renormalization group irreversibility. One of the implications of 
this proposal is 

Sain) > S a (r 2 ) for n > r 2 . (54) 

This means c a {r) > in our notation. We note that, at least for the entanglement entropy, 
this relation follows from translation invariance only while (|53|) with a = 1 requires 
the full Poincare group symmetry, which is an essential ingredient for the c-theorem |23j . 

6 Appendix A: Correlator of exponential operators in the 
sine- Gordon model 

The two-point correlation function of (|3()j) can be expanded as a sum over form factors as 
follows [15] 

/. e is/4^a<t>(r) ,. e ivW<6(0) A = y J_ du 1 ...du 2n [ TT e _2 ^ v" i+ "i) ] 

x 1 jsfao J ° Vf=i / 

Xf a (ut, ...,U2n)/a'( n 2n, , (55) 

where f a (ui, ...,u 2n ) is the form factor 

/.(til, ...,«*») = (-l)^)/ 2 (^|^) n (ft (^)") >< ...,«*»), (56) 
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and 

A{u 1 ,...,U2n) = J ~ =^ , , • (57) 



Ut<j <n ( U i ~ U i) Un+l<i<j ( U i ~ U j) 
Ur=l Ulln+l i u r + ««) 



This series can be expressed as a Fredholm determinant |15| 

'. e iV^a<t>(r) .. e iv^a'0(O) = _ \ 2 R a ^ a , R^_ a ,) , (58) 

where Rg is an integral operator on the half-line (0, oo) with kernel 

(-) " T . (59) 

Rj is the transpose of Rg, x = mr, and 

A = — (sin(7ra) sin(7ra')) 5 . (60) 

Defining 

r = log det(l - X 2 RgRj) (61) 
we have the following differential equations [2U] 

d 2 r 1 dr , 9 , 

ax z x ox 

with boundary condition 

d(x, A) ~ 2Ai^g(x) as x — > oo . (64) 

In the case relevant for this work a' = —a, 9 = 2a and 

i 

A = — sin(-7ra) (65) 

7T 

is purely imaginary. This translates into a purely imaginary d. Thus, defining v a = —id 
equations ifST]). (Sg) and (J33J) follow. 

Equation (j63j) has appeared in the literature in different forms. It adopts the standard 
Painleve V form [21] by the transformation 

i = ii±£, (66) 

which leads to 

" 1,' /l 1 \ / '\ 2 e 2 (h-l) 2 / 1\ + . . 

+ -/i = — + J ( /i ) + — =-£- U - - J + 2 /, / . 67 

t V2/i h-lj\ J 8 x 2 V V (^-!) 

Also, by the transformation = sinh(/) we get the equation 

," 1 ,' 1 w « # 2 tanh(/) , . 
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which differs however from the equation given in |15j by a factor 4 in the last term. It can 
be checked that the well-known differential equations for the Ising model spin and disor- 
der correlators [211 125j . (cr(r)er(O)) and (p(r)p(O)), can be obtained from the differential 
equations (|62j) and (|63j) for the sine-Gordon correlators through the identification 26 



(a(r)a(O)) 2 = ^sin(^ y/^<j>{r)) sin(^ y/^<f>(0))j , (69) 
(p(r)p(O)) 2 = ^cos(i\/47r^(r))cos(i\/47r^(0))^ . (70) 

7 Appendix B: Numerical evaluation in the lattice 

The vacuum expectation value {Oa} for any operator Oa localized inside a region A must 
coincide with tr(p^O^), where pais the local density matrix. This fact was used in [2Z] 
to give an expression for pa in terms of correlators for free Boson and Fermion discrete 
systems. We use this method here to compute the entanglement entropy for a free Dirac 
field on a lattice. 

Consider a lattice Hamiltonian of the form 

H = Y f M ij 4c j , (71) 

where the creation and annihilation operators c|, c, satisfy the anticommutation relations 

{c\,Cj} = Sij. Let the correlator Cjj = (cjcj^ be taken on any given fixed eigenvector of 

the Hamiltonian and call C{\ to 
density matrix pa has the form 



the Hamiltonian and call C„4 to the correlator matrix restricted to A . Then the reduced 



where the di are independent fermion annihilation operators which can be expressed as 
linear combination of the C{ and cj. The e/ are related to the eigenvalues v\ of C A by 

"' (73) 



\-v l 



Hence both the entropy and logtrp n can be evaluated as Sa = S S; and logtr// 1 = £ 
logtrp", where 

S t = log(i + e~ tl ) + e, e /\ = -(1 - ^)log(l - n) ~ \og(n)n , (74) 

1 + e ei 

logirpf = log(l + e- ne 0-^log(l + e" ei ) = log((l-^) n + ^ n ). (75) 

We are interested in the vacuum (half filled) state of a Hamiltonian with symmetric spec- 
trum around the origin. In this case 

C = 9(-M) , (76) 

where 6(x) = (l+sign(x))/2. 
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The lattice Hamiltonian for a Dirac fermion follows from the discretization of the 
Hamiltonian for a Majorana field, to cope with fermion doubling. This is given by 

. N-1 N-1 
W=yE + *n*n+l) « , (77) 

n=0 n=0 

with tfj, = ^ = ¥* and {* ^ = 

fimnfiij ■ We have taken N lattice sites and set 
the lattice spacing to one. Redefining the fermionic operators as 



C2k 



l(*lk+Mlk), (78) 



with {cn, c m } = 5 mn we have a Hamiltonian in the form (|71|) 

. AT-1 JV-1 
2~ 



n = T E ( C n+l C « " C ^n+l) + m £ (-l)Mc„ (80) 



n=0 n=0 

with 



Mij = —(5(i,j - 1) - j + 1)) +mi(ij)(-l) n . (81) 



-2 

2~ 

The correlators for the infinite lattice limit N — > oo are then 



/ t \ If / /" 2 r "T.cos(27rx(i — j)) . . . ,„_,. 

<4 c i> = o%-i),o + (-l) / dx 7 r F== # for z-j even, (82) 
z jo ym^ + sm(27rxj z 

(etc,-) = === = sin(27rx(i - j)) fori-jodd. (83) 
jo v »n + sm(27rx) / 

Similar expressions where used in [8]. 

From (|74|) and (|75|) we can obtain c n (r) and c(r) by taking numerical derivatives. We 
have used 

/' (<1 + I) = 7 (/(« + 2 ) + /(« + 1) - /(?) " /(? " !)) • ( 84 ) 



This form for the derivative smoothes small oscillations which appear for quantities de- 
pending on a lattice size q when evaluated between even an odd adjacent values of q. 
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